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Abstract 

^ I The AMS-02 collaboration has recently released data on the positron fraction e''"/(e + 

2^ I e"*") up to energies of about 350 GeV. If one insists on interpreting the observed excess as 

fT^ ' a dark inatter signal, then we find it is best described by a TeV-scale dark matter anni- 

hilating into T^T~ , although this situation is already severely constrained by gamma-ray 
^T . measurements. The annihilation into iJi^ lJi~ is allowed by gamma-rays more than t~^t~ ^ 

^f^ • but it gives a poorer fit to AMS-02 data. Moreover, since electroweak corrections in- 

duce correlations among the fluxes of stable particles from dark matter annihilations, the 
recent AMS-02 data imply a well-defined prediction for the correlated flux of antipro- 
^ I tons. Under the assumption that their future measurements will not show any antiproton 

excess above the background, the dark matter interpretation of the positron rise will pos- 
sibly be ruled out by only making use of data from a single experiment. This work is 
the first of a program where we emphasize the role of correlations among dark matter 
signals. 
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1 Introduction 

Despite the overwhelming evidence that Dark Matter (DM) contributes to a significant fraction of 
the energy density of the universe, its nature and properties are stih unknown. Many experimental 
activities are nowadays performed assuming the particle nature of DM. In particular, in the indirect 
searches of DM [1] the goal is to detect the products of DM annihilations or decays around the Milky 
Way: Standard Model (SM) particles are produced as primaries and they subsequently decay into SM 
stable states such as electrons and positrons, protons and antiprotons, deuterium and antideuterium, 
photons (7-rays, X-rays, synchrotron radiation) and neutrinos. In praticular, the recent positron 
data released by AMS-02 [2] may be interpreted as due to DM annihilations. 

The necessary inclusion of ElectroWeak (EW) corrections [3, 4, 5, 6] significantly alters the final 
spectra when the mass of the DM particles is larger than the EW scale (see also [7, 8, 9]): soft EW 
gauge bosons are copiously radiated, thus opening new channels in the final states which otherwise 
would be forbidden if such corrections were neglected. As a consequence, all stable particles are 
present in the final spectrum, independently of the primary channel of DM annihilation/decay. 
An inevitable consequence of this is inducing correlations among the predicted fluxes of SM stable 
particles targeted by the various experiments. For example, within a given model of DM, EW 
corrections to DM annihilations/decays lead to correlated fluxes of antiprotons (p) and neutrinos 
whose detection is one of the goals of AMS-02 and IceCube, respectively. 

This work is the first of a program where we investigate the role of correlations among DM signals. 
Based on the idea that EW interactions connect all particles of the SM and that correlations can be 
established among different fluxes, our logic is very simple: 

1. suppose a signal is detected in a particular channel and assume it is of DM origin; 

2. then the EW corrections imply that the fluxes of other cosmic-ray species should be modified 
accordingly, in a very specific way depending on the nature of the process which is assumed in 
part 1; 

3. these correlations can then be tested against existing data, or used as predictions for upcoming 
data, to verify or reject the initial hypothesis of DM interpretation of the signal. 

This approach is essential to draw robust conclusions about DM origin of a signal, and it is particu- 
larly relevant now, due to the amount of upcoming data from very different sides and channels (e.g. 
AMS-02, Fermi, HESS-2, IceCube, Pamela). 

In this paper we apply the general logic described above to the specific case of the AMS-02 
experiment, which has recently released positron data [2] and will release data from anti-proton 
measurements in the very near future. We first consider the recent AMS-02 data on the positron 
fraction e^ ji^e^ + e~) (which extend to higher energies the observations by Pamela [10] and Fermi 
[11]) and fit them as a signal of DM. Then, the EW corrections provide a prediction for the correlated 
flux of p (which are already constrained by current Pamela data [12]). Once this prediction is made, 
we may compare it with the projected sensitivity of AMS-02 for the p channel under the assumption 
that the measurements of AMS-02 do not show a signiflcant p excess above the background. 

The correlation between e^ and p is particularly relevant as it is a very clear example to illustrate 
the role played by EW corrections in constraining DM models even within a single experiment, such 
as AMS-02, which is able to measure the fluxes of different particle species. The inclusion of EW 



corrections will allow us to answer a simple and interesting question: will it be possible to rule out 
conclusively the DM interpretation of the positron rise using only the data collected by the AMS- 
02 experiment (and thus disregarding other experiments)? Notice that the correlations between 
different final states can also help to break the degeneracy implicit in the interpretation of the e^ 
signal from either DM or astrophysical sources, mainly pulsars. While it was recently concluded 
that future e data will likely be insufficient to discriminate between the DM and the single pulsar 
interpretations of the cosmic-ray lepton excess [13], this conclusion did not take into account EW 
corrections and the correlations induced by them. The correspondingly generated flux of p might be 
detected by AMS-02 and not be ascribable to pulsars. We will come back to this point in a separate 
publication [14]. 

The paper is organized as follows. In Section 2 we briefly discuss the signals and the backgrounds 
relevant for our analysis. Our results are presented and discussed in Section 3, including a description 
of the salient features of our detector simulation. The concluding remarks are collected in Section 4. 

2 Signals and backgrounds 

As a signal, accounting for the positron excess observed in the data, we consider DM annihilations 
in the following leptonic channels 

DM DM ^e+e", ^+/x~, r+r", (2.1) 

each with 100% branching ratio. We will not consider the case of decaying DM. The above anni- 
hilation channels are the best situations for interpreting the positron rise as a signal of DM. Other 
non-leptonic annihilation channels, such as bb and VF"*"!^", easily produce hadrons in the final state, 
e.g. p, and are already excluded by the p Pamela data [12], using EW corrections [3]. Note however 
that other complimentary data put already strong constraints on these "leptophilic" scenarios, e.g. 
7-rays [15, 16]; as mentioned above, we insist on exploring the exclusion capabilities of data from a 
single experiment, with EW corrections taken into account. 

The correlations we want to study in this paper are between positron fraction and p data, 
therefore the Cosmic-Ray (CR) species we have to deal with are e~ , e^ ,p. An example of how DM 
can account for the positron excess without upsetting the p flux is shown in Fig. 1. In the following 
we will only consider the data points with energies above 15 GeV, in order to be totally sure that the 
fluxes are not affected by solar modulation. To generate the signal fluxes at production, i.e. positrons 
and p from DM annihilations, we use the numerical code PPPC4DMID described in Refs. [17, 3], 
which includes EW bremsstrahlung from external legs in the factorization limit. The predictions of 
cosmic-ray fluxes originated from DM annihilations crucially depends on how the DM is distributed 
in the galactic halo. This astrophysical uncertainty is currently irreducible, therefore we consider 
three of the usual profile choices: 
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(2.2) 
As for backgrounds, "primary" electrons can come from galactic CR while interactions of CRs with 
the interstellar medium is a source of "secondary" electrons, positrons and p. The propagation of 
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Figure 1: Dark matter interpretation of the excess in AMS-02 positron fraction data (left panel); the 
corresponding flux of p would not show significant deviations in AMS-02 p mock data (right panel), generated 
according what discussed in Sect. 3.2. The dashed blue lines correspond to the contribution of the astrophysical 
backgrounds, computed with method 2 described in Sect. 2.2. The solid black line is the total flux including 
the DM signal, where we have chosen DM DM ^- t^t~ , Mdm = 1 TeV, av — 2.5 x lO^^'^cm'^s"^, Burkert 
halo profile. 

the signal and background fluxes from their production region to the detector is affected by many 
processes, mainly diffusion and energy losses. The number density per unit energy nj(r, z,p) of the 
cosmic ray species of interest i (= e"'",e~,p), at any given point of cylindrical coordinates r,z and 
momentum p, evolves according to the transport equation [21] 
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where D {Dpp) is the diffusion coefficient in position (momentum) space, Q{r, z,p) is the source term 
due to DM annihilations, which depends on the DM mass, annihilation cross section and the energy 
spectrum at the interaction point (before propagation), Vc is the convection velocity, Tsp is the time 
scale for nuclear spallation processes, tj is the time scale for fragmentation. The flux is then given 
by the solution of the transport equation through $j = /3jni/(47r), where /3j is the particle velocity 
in units of the speed of light. 

The coefficient Dpp determines the diffusive reacceleration and is proportional to the the inverse 
of D 

D -1 ^ 1 (2 4) 

where va is the Alfven velocity. Eq. (2.3) is solved in steady state [dn/dt = 0), in a cylinder of 
height 2L and radius i? = 20 kpc, encompassing the galactic plane. The diffusion coefficient is 
conventionally parameterized as D{E) = I3Dq{TZ/1GV) , where TZ its rigidity. Both the position 
dependence and the deviation from the linear dependence on /3 have been neglected. For the profile 
of the galactic magnetic field, which is rather uncertain, we choose [22] 



B{r,z) = Bqc 10 kpc 2 kpc 



(2.5) 



where r© = 8.3 kpc and Bq is of the order of (1 — 10) /iG. 

We carry out our analysis following two different methods for handling the propagation of signal 
and background fluxes: in method 1, we take a conservative approach letting the astrophysical 
uncertainties largely vary; in method 2, we assume a specific propagation model and use it to 
consistently propagate the signal and background fluxes of all relevant species. Strictly speaking, 
the correct way to proceed would be to choose a specific CR propagation model within which to 
compute both signal and background for all particle species and then scan over the model parameters. 
We decided to adopt these two different methods because of simplicity and because versions of each of 
them are often used in the literature. They are meant to represent two extreme ways of proceeding, 
so the "truth" is very likely to lie in between of the results obtained with these two methods. 

2.1 Method 1 

The signal fluxes of e and p due to DM annihilations, propagated to Earth, are computed with the 
code PPPC4DMID, for different halo profiles and using a semi-analytic propagation model known as 
"MED" [23, 24]. This model corresponds to a specific choice of the coefficients of the diffusion-loss 
equation, for both e± and p: L = A kpc. Do = 3.38 x lO^'^cmVs, S = 0.70, T4 = 12km/s, Va = 52.9 
km/s. Bo = 4.78/iG. 

As a reference background p we use a spectrum which is between the Kolmogorov (KOL) and 
Kraichnan (KRA) models of Ref. [25] (the same choice has also been done in Ref. [26]). Changing 
the propagation model does not significantly affect the shape of the p flux, but it changes the 
normalization [27]. For primary electrons and secondary e , we use the parametrization [28], also 
used in [29] (see instead Ref. [30] for a more data-driven approach) 
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where e = E/ GeV. Each of these background functions suffer from astrophysical uncertainties 
related to the CR propagation in our Galaxy. We deal with theses uncertainties by simply allowing 
the backgrounds to have floating normalizations and slopes; in practice, we multiply each of the 
background fluxes by Ai e^' (z = e~ ,e'^,p), and marginalizing over Ai £ [0.5, 2] and pi £ [—0.05, 0.05]. 
This approach has the virtue of accounting for the astrophysical uncertainties of the propagation 
model in a simple and conservative way. It will therefore result in rather cautious bounds on the 
parameter space. On the other hand, the relationship between the A,p and the physical parameters 
of the transport equation is obscured. Furthermore, the marginalization over A,p parameters is 
done independently for e and p. This corresponds to treat the normalizations and slopes of the 
backgrounds of these types of fluxes as uncorrelated, as they are propagated with different physical 
parameters. 

2.2 Method 2 

Alternatively, we consider another way to deal with astrophysics. We insist on the fact that we 
want to consider a consistent setup where both signal and background are propagated with the same 
set of propagation parameters. So we consider a specific propagation model defined by L = 4 kpc. 
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Figure 2: Fits of our reference propagation model, as defined in Sect. 2.2, to current data: Pamela p (left 
panel), Pamela p/p (right panel). Dot-dashed lines refers to the case without correcting for solar modulation 
effects. 

Do = 1.50 X lO^^cmVs, 6 = 0.50, Va = 13.715 km/s, Bq = 7.0/iG. Unlike the MIN, MED, MAX 
models, no convection is considered. We set the injection index of primary protons to be 2.35, while 
the injection index of primary electron spectra is 0.6 (2.53) below (above) the reference rigidity of 
4 GV. This model provides a good fit to[ll], Pamela positron [10] and p [12] data (see Fig. 2), as 
well as to Fermi all-electron data at low energy and B/C data. 

The signal fluxes due to DM annihilations are computed by computing the energy spectra at 
the production point with the PPPC4DMID code, and then propagating them using GALPROP 
[31] and the model specified above. The background fluxes for e^ and p are also computed using 
GALPROP, for the same propagation model used for the signal. 

Notice that, differently from method 1, for method 2 we assume we have a rather precise knowl- 
edge of the astrophysical details. As of today, this is not the case, but for instance the future data 
releases of AMS-02 on absolute e~^,e~ ,p,p as well as light nuclei ratios, will improve the accuracy 
of the propagation parameters and of the resulting backgrounds (see e.g. Refs. [24, 32]). Therefore, 
carrying on this method serves as a representative way for what can be achieved with more detailed 
information about CR propagation in our Galaxy. 

This approach has the virtue of propagating consistently all fluxes within the same model, al- 
though it pays the price of not being generic, as a scan over the uncertainties of the model parameters 
is not performed. 



3 Results 

3.1 Analysis of the positron fraction data 

We flrst perform a x^ analysis of the AMS-02 positron fraction data for the different channels under 
consideration e+e~, //+/i~, T"'"r~, in order to identify the best-fit one. We also checked that best-fit 
regions obtained by fitting the positron signal with the bb and W~^W~ channels are already exlcuded 
by Pamela anti-p data [12]. As already mentioned, we restrict ourselves to use only the E > 15 
GeV bins of the data, as the lower energy bins are strongly affected by solar modulation effects, so 
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Figure 3: The x function as a function of the DM mass, of the AMS-02 positron fraction data, for method 
1 (left panel) and 2 (right panel). The halo profile is Einasto. Different colours correspond to different DM 
annihilation channels: /^^/i~ (blue), t^t^ (red). The e'^e' channel gives an even higher x^ , f^ot shown in 
the plot. 

a total of 36 data points. The numbers of degrees of freedom are 36-6=30 and 36-2=34, for method 
1 and 2, respectively. 

In Fig. 3 we show the x^ as a function of the mass, minimized over the annihilation cross section, 
for the /i and r channels. As a reference halo profile we chose Einasto, but for other choices the 
results are very similar. The e channel is not shown as it gives a very poor fit with even higher values 
of x^: x^/dof > 6.6, 5.8, for methods 1 and 2, respectively. The /i channel gives Xmin/^of — 1-9, 2.4 
for methods 1 and 2, respectively, so it also provides a rather poor fit to the data. Instead, the r 
channel gives Xmin/^of — 0.7, 1.0 for methods 1 and 2, respectively. So we conclude that the only 
case producing a good fit to the AMS-02 positron fraction data are DM of about 1 TeV annihilating 
into t'^t" . We also note that the case of DM annihilation into /i"*"//", although it gives higher x^ 
to AMS-02 positron data than t~^t~ ^ it is less constrained by gamma-ray measurements [33]. 

After having identified the situation which fits the signal at best, we now want to proceed and 
correlate it with future measurements of anti-p flux by the same experiment. To this end, we simulate 
the detector response and generate mock data. 

3.2 Simulation of AMS-02 p data 

As we have mentioned in the Introduction, we work under the hypothesis that the future data on 
the p flux reported by AMS-02 will not show any evidence of DM annihilations and will be therefore 
compatible with the expected background. In order to assess the projected sensitivity of AMS-02, 
we assume that the p flux will be measured up to energies of the order of 300 GeV and that the 
number of detected p in a, given kinetic energy bin (.Ekin — A£^kin/2, -Ekin + A£^kin/2) is 

A^p-(^kin) - ap(^km)^p-(^kin)ASkmAt , (3.1) 

being ap the geometrical acceptance, ^p the p flux and At the exposure time, and we have taken 100% 
efficiency [34]. Assuming a Poisson distribution, the statistical error is therefore (AA'p)stat/-^p = 
1/ yJNp. For the systematic error for p detection and reconstruction, we can only suppose the syst 
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Figure 4: The red contours are the best-fit 3a regions of the parameter space [m^, (o'^')); of the AMS-02 
positron fraction data with DIM DM — > t^t^ channel, for method 1 (left column) and method 2 (right 
column), as described in the text. The blue lines are the projected exclusions (at 99.7% CL) due to p AMS- 
02 mock data, for 1-year (thin blue line) and 3-year of data taking (thick blue line). From top to bottom: 
Burkert, Einasto, NFW profiles. 



error is in the 1-10% range [35]; therefore we consider as reasonable values {ANp)syst/Np = 5% 
for the 1-year data and {ANp)syst/]^p = 3% for the 3-year data, to be added in quadrature to the 



statistical error. For the geometric acceptance of p we choose ap = 0.2 m^ sr as a safe lower limit 
above 1 GeV [36]. Finally, we use a linear approximation of the detector energy resolution (see 
[37, 26]), A^kin/^kin = (0.042 • (Ekin/GeV) + 10)%, which allows about 20 bins per decade. 

Within this setup, we generated mock p data for AMS-02 in the range (10 — 300) GeV of kinetic 
energies, following the two reference backgrounds discussed in the previous section. In the right panel 
of Fig. 1 we show the mock data generated as described above, as well as the reference background 
and the DM contribution to the p flux, for method 2. 

3.3 Correlations between e"*" and p data 

For the best-fit situation DM DM — ?• t^t" we find the 3a best-fit contours, on the (ttidMi (o"^)) 
plane, corresponding to Ax^ < 11.8, of the AMS-02 positron fraction data. They are shown as 
red contours in Fig. 4, for the two methods discussed in the previous section. On the same plots, 
the blue lines are the projected exclusions (at 99.7% CL) due to p AMS-02 mock data, for 1-year 
and 3-year of data taking, generated according what discussed in the subsection 3.2, and assuming 
the data will not show deviations from the background. We checked that varying the geometrical 
acceptance ap by a factor of 2 does not change the results significantly, as the errors are dominated 
by systematics. We vary the DM profile, following the choices in Eq. (2.2), and we use the two 
different methods described in Section 2. We reiterate that significant portions of the parameter 
space shown in Fig. 4 are already excluded by complementary data, such as 7-rays [33]. However, 
the bounds for shallow halo profiles like isothermal or Burkert are weaker. We do not include these 
constraints in our plots as the focus of this paper is on the correlations between AMS-02 current 
and projected data. 

The results are not altered significantly by choosing Einasto or NFW halo profiles, but exclusions 
are weaker for the Burkert profile, as it assumes lower DM density near the center, and therefore a 
lower p signal would be produced. 

Furthermore, the results are sensitive to the method chosen for handling the propagation of 
signal and background. The fact that the best-fit contours of the positron fraction data are smaller 
for method 2 than for method 1 is not a surprise, as in method 2 we have less freedom to vary the 
astrophysical background and therefore the DM signal is more constrained. For the same reason, the 
exclusion curves coming from p projected data are more stringent for method 2. Therefore, adopting 
method 1, any conclusion is weaker than for method 2, and we can imagine the actual situation will 
lie in between of these two scenarios. 

We can conclude that 1-year of AMS-02 data on p would already be enough to almost rule 
out the DM interpretation of the positron excess. With 3 years of data all the exclusion curves are 
obviously more stringent. 

The high energy bins of the AMS-02 positron fraction data show a decrease in the slope and 
suggest a possible turning down of the signal. If future positron data releases will show more clearly 
a fall of the positron fraction, the best-fit contours with DM signal would shrink and better pin 
down the DM mass. 



4 Conclusions 

This paper demonstrates the relevance of correlations induced by the EW corrections in indirect DM 
searches. The results of the analysis described above allow us to draw the following conclusions: 

• if the recently released AMS-02 positron fraction data are interpreted as a signal of annihilating 
DM, then the x^ analysis favours a a DM of about 1 TeV annihilating into t'^t^ (a situation 
which is already severely constrained by 7-ray measurements); 

• assuming no signal is seen in future p data released by AMS-02, one will be able to exclude 
almost completely the DM interpretation of the positron excess, by using only data from a 
single experiment. 

These conclusions have to be considered together with two caveats. First, while the energy spectra 
of SM particles at the interaction point coming from DM are predictable in terms of particle physics, 
the propagation of these particles in the complex galactic environment is still uncertain. Therefore, 
astrophysical uncertainties and the way to deal with them affect the results. We have pointed out 
the importance to use consistently the same propagation setup for both signal and background, 
which is especially relevant when combining together the fluxes of more than one particle species. 
Second, although the assumptions we made are consistent with the understanding of the AMS-02 
detector from outside the collaboration, our results may be sensitive to the detector features we used 
for generating the projected data. 

This work should be regarded as the first of a program where we will investigate more aspects 
and examples of correlations among DM signals. In subsequent publications we will study different 
channels and update these results as soon as new data become available. 
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